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ABSTRACT 


The  solutions  of  the  phase  equation  for  a  potential  with  the 
asymptotic  properties  of  an  atomic  polarization  potential  are 
studied  with  the  intention  of  developing  a  computer  program  for 
direct  integration  of  the  phase  equation  for  the  s-wave  phase 
shift  of  an  arbitrary  potential.  Such  a  program  is  developed  and 
its  limitations  are  discussed.  In  addition,  an  appendix  is  devoted 
to  the  variational  principles  of  Kohn  and  Schwinger  for  calculating 
phase  shifts. 


LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL  DUDLEY  KNOX  LIBRARY 

MONTEREY.    CALIF.    93940  NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA  93943-5101 

TABLE  OF  CONTENTS 

I.      INTRODUCTION 9 

II.      THEORY i:L 

A.  THE  PHASE  EQUATION n 

B.  THE  POTENTIAL 13 

III.   CALCULATIONS 16 

IV.   RESULTS 22 

V.   CONCLUSIONS 28 

APPENDIX  A.l  -  THE  KOHN-HULTHEN  METHOD 29 

APPENDIX  A.  2  -  THE  SCHWINGER  METHOD 33 

COMPUTER  PROGRAM 38 

REFERENCES 44 

INITIAL  DISTRIBUTION  LIST 46 

FORM  DD  1473 47 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Plot  of  the  Potential  Function 15 

2.  Plot  of  Phase  Function  for  K  =  0.1,  Yukawa  Potential  18 

3.  Plot  of  Phase  Function  for  K  =  1.0,  Yukawa  Potential  19 

4.  Plot  of  Phase  Function  for  K  =  0.1,  Yukawa  Potential 20 

5.  Plot  of  Phase  Function  for  K  =  0.1,  Polarization 
Potential 23 

6.  Plot  of  Phase  Function  for  K  =  1.0,  Polarization 
Potential 24 

7.  Plot  of  Phase  Shift  vs.  K  for  Polarization 

Potential 27 


SYMBOLS 

r  THE  RADIAL  COORDINATE  IN  BOHR  RADII 

K  THE  ENERGY  IN  RYDBERGS 

n  (r,K)  THE  PHASE-FUNCTION  IN  RADIANS 
o 

n  (K)  THE  PHASE  SHIFT,  n  (°°,K)  ,  IN  RADIANS 

u(r)  THE  SCHRODINGER  RADIAL  WAVE  FUNCTION 

U(r)  orV(r)  THE  POTENTIAL  FUNCTION  IN  RYDBERGS 


ACKNOWLEDGEMENTS 

The  writer  wishes  to  express  his  appreciation  to  Professor 
R.  A.  Armstead  for  his  understanding,  encouragement,  and  aid 
during  the  course  of  this  investigation  and  also  to  my  wife  Pauline 
for  her  unquestioned  support,  loyality,  and  gentle  goading  as  well 
as  her  skill  with  the  typewriter. 


I.   INTRODUCTION 

The  Quantum  Theory  of  scattering  processes  is  of  basic  import- 
ance in  studies  of  atomic,  nuclear,  and  elementary  particle  physics. 
Analysis  of  Atomic  collision  processes  is  required  in  such  diverse 
applications  as  gas  dynamics,  chemical  reaction  theory,  and  stellar 
dynamics.  Nuclear  Physics  relies  on  scattering  of  various  particles 
from  nuclei  for  a  great  deal  of  the  quantitative  information  about 
nuclear  processes  while  in  Elementary  Particle  physics,  very  little 
other  information  is  available  except  for  the  various  scattering 
events  that  are  observed. 

There  are  many  approaches  to  the  scattering  problem,  but  one  of 
the  most  useful  (especially  in  low  energy  applications)  is  the  method 
of  partial  wave  analysis.  The  method  was  first  introduced  by  Faxen 
and  Holtzmark  based  on  a  method  developed  by  Lord  Rayleigh.  Rather 
then  delve  into  the  basis  of  the  theory,  it  will  be  assumed  that  the 
reader  is  familiar  with  the  treatment  of  Schiff  [Ref .  1]  or  Merzbacher 
[Ref .  2] . 

A  particularly  interesting  procedure  for  calculating  scattering 
phase  shifts  is  the  "Phase  Method."  Originally  introduced  by  Courant 
and  Hilbert  [Ref.  3]  and  applied  to  scattering  by  Morse  and  Allis 
[Ref.  4] ,  the  method  has  recently  been  revived  and  expanded  in  a 
series  of  papers  by  Calogero  [Ref.  5,  6]  and  his  new  book  [Ref.  7] . 

This  paper  will  investigate  the  application  of  the  phase  equation 

to  a  potential  which  has  the  asymptotic  properties  of  an  atomic 

4 
polarization  potential  (1/r  ) .  The  method  employed  is  a  direct 


integration  of  phase  equation  (11) .  In  the  appendix,  the  variational 
principles  of  Kohn  and  Schwinger  are  treated  and  are  suggested  as 
alternate  methods  for  computation  of  phase  shifts. 
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II.   THEORY 

A.   THE  PHASE  FUNCTION  METHOD 

The  phase  function  method  is  a  method  of  reducing  the  second  order 
linear  Schrodinger  Equation  to  a  first  order,  non-linear  Riccati  equa- 
tion. For  s-wave  scattering,  the  solution  for  phase  shifts  reduces 
to  the  solution  of  the  differential  equation 

(1)        ^(r)     =   -  [V(r)/K\   SiH20<r  t  y\) 

To  show  how  this  equation  comes  about,  it  is  convenient  to  approach 
the  problem  from  the  standpoint  of  A.  Ronveaux  [Ref .  8] .  He  starts  by 
considering  a  second-order  homogeneous  differential  equation  with  two 
perturbative  terms ,  V,  (X)  and  V-  (X)  ; 

«    u"  -t    flt*W  *   Q(x)u   -  V,(y)u  f  Vzu)  u' 

with  u(0)  =  y  and  u'  (0)  =  v.  This  can  be  written  in  the  form  of  the 
following  system 

<3)        X'(*>   r  M*>  X(x)    t   v(x)  x(x) 

where  /        K  . 

X(x)    -  (u\  f(o]   = 

A(x)  -■/  o       i  \        V(x)  ■-  (o     O 

(-Q  -?)  (  v,  v2 

The  unperturbed  solution,  X  is, 

(4)         x  (*)   -    w(*>  C 
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a,      u2 


where     C   -    /  c 

\C2. 
is  a  constant  vector,  and 

W(x)    - 

U,      U; 

and  u,  and  u~  are  two  linearly  independent  solutions.  The  solution  of 
the  complete  system  may  then  be  written 

(5)         t(x)  ■■     W(x)CCx-) 

Using  the  Lagrange  Method  of  Variation  of  Constants,  then 

C'OO    -     j  W(x)J"lVM  VJ(x)  C(*) 

or 

C'(k)    r    rA(x)C(x) 

where 


C(o)  --   ^  WCo)|  /^\ 


If  we  now  define  the  phase  function 

W  SU)  -  C4Cx)/C,(<) 

and  differentiate,  we  find  that 

(7>     s'(x)  --  M7iU)  +[rtZ2(x)-rill(x)]s(x)  -fin(xys2M 


Evaluating  M  gives  the  phase  equation 
(8) 


S'(x)   =    |[w,U)]    *  f  u2u)  S(x)]j  /dt-i  WU) 


where  det(W)^  0  and 


[J    u>)/<j+  [<V«)vj[ 
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Given  the  general  form  one  now  approaches  the  problem  of  s-wave 
phase  shifts.  The  Schrodinger  Equation  (in  natural  units,  n  -C  -  I  ) 
is 
(9)     Ui'Xt)      +  Kl  u(r)     *      V(r)  u(r) 

choose 

UCo^    ~  /*.  -    O 

Then  the  general  solution  is 

U(r)    =    C,     S/n  kr    -  C2   cos  l<r 

r^       C'fCr)    [sinKr    -    SCr)    Co^^] 
or 

UCr)    ~       C,  (r)  cot  rf0(rtK)   Sm  (  Kr  +  ty) 

where 

S-Crr)  =    c2<r)/c,rr;=     -  "fan  ^(f,)c) 
Then,  letting   V,  -  V,  V2  =  0}U,-  s/nkr}  uz~  c<?s/<r,  cMk/-   *     the  general 
phase  equation  becomes 


(10) 


7 

(d/dr  )  +an  YJ0      =      -    [  V(r)  /k}(  £•**!  kr    f    fa  n  ^/os  /<rr  J 


with    +aY>Yfe  (o,l()zo.      Noting  that      (  d  ld<  )  ion  V)0    :    S<ec    ^(d  *]<,/df  ) 
it  becomes  apparent  that 


(11) 


(d/dr)T[fl    =       -  [v(r)/l<]  Sin2  iKr  t  vji«  ) 


B.   THE  POTENTIAL  FUNCTION 

The  object  of  this  paper  is  to  investigate  the  properties  of  the 
numerical  solutions  for  the  Phase  Equation  (11)  when  applied  to  a 
potential  which  has  the  asymptotic  properties  at  infinity  of  an 
Atomic  Polarization  potential ,  i.e., 

U(r)    >  (const)  x    '  M 

r — 7  &* 
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In  order  to  simplify  the  analysis,  an  additional  condition  was  imposed 
on  the  potential  in  the  form  of  a  requirement  that  as  r_>  O  ,  the 
potential  should  approach  zero  as  if  ,  i.e., 

U(f)    __>  (con&i  )xr 

V — ••?  O 

The  function  chosen  to  satisfy  these  conditions  was 


The  asymptotic  properties  of  this  form  are  obvious,  however,  it  is 
difficult  to  calculate  with.  In  order  to  simplify  the  problem  further, 
the  substitution  T-  xr       was  chosen  which  leaves 

(12)  U(r)    -~      r  idU  xv  e 

Jo 

It  is  no  longer  readily  apparent  that  U(r)  has  the  correct  asymptotic 

properties,  but  it  may  be  easily  verified  by  integration  and  power  series 
expansion.  The  success  of  this  form  is  shown  in  Fig.  1. 
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V(r)   -r 


0.08  — 


0.04 


Figure  1 
Plot  of  the  Potential  Function 
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III.   CALCULATIONS 

Having  derived  the  phase  equation  (11) ,  one  now  wishes  to  devise 
a  method  to  put  it  to  practical  use.  Calogero  develops  a  variational 
scheme  for  the  phase  equation  [Ref .  9]  while  other  authors  deal  with 
the  equation  in  other  ways  [Ref.  10,  11]  .  An  excellent  review  has 
been  prepared  by  Babikov  [Ref.  12]  which  provides  a  good  general 
introduction  to  actual  calculations  and  an  equally  good  bibliography. 
Cwing  to  the  fact  that  the  equation  is  non-linear,  the  various  ap- 
proaches have  dealt  with  some  form  of  approximate  solution.  This 
paper  will  deal  with  the  relatively  simple  problem  of  constructing 
a  program  to  directly  integrate  the  phase  equation  using  numerical 
approximation  only. 

The  particular  method  chosen  was  a  fourth  order  Runge-Kutta 
Integration  scheme  to  be  implemented  on  the  IBM- 360  time  sharing 
system.  General  discussions  of  the  numerical  method  may  be  found 
in  numerous  places  [Ref.  13,  14] ,  and  so  only  the  necessary  portion 
of  that  theory  is  presented  here. 

The  Runge-Kutta  method  proposes  to  solve  the  equation 

(i3)  y/  -     {(*>  y) 

subject  to  the  initial  conditions  X  =  X  when  Y  =  Y  .  The  iterative 
J  o         o 

solution  process  proceeds  via  the  following  equation; 

Vri+,     -     y,    +    1  I  k,  4  2K2    +  ZKz    +   *m] 

(d 
16 


(14) 

where 


|<4  -       l\  $  (  X,  *  h  >  y>  +  ft  ) 


and  where  X0+ 1  =  *>,  t  )i   .  The  programming  needed  to  provide  a 
solution  of  the  equation  (11)  which  gives  the  exact  phase  shift, 
yj0(.K)  -  V0(^°)  k)  /  then  reduces  to  programming  a  simple  iterative 
procedure.  The  program  developed  is  printed  at  the  back  and  is  in 
such  a  form  that  it  may  be  run  either  on  the  time-sharing  system  for 
quick  approximations  or  directly  on  the  computer  system  in  the  regular 
batch  processing.  A  second  version  of  the  program  used  to  obtain  a 
print  out  of  the  values  of  the  phase  function  at  various  r  values  is 
also  included.  The  use  of  the  second  form  will  become  obvious  as 
the  discussion  unfolds. 

It  now  becomes  necessary  to  consider  the  actual  running  of  the 
program  and  its  particular  limitations.  In  order  to  properly  investi- 
gate the  program,  the  Yukawa  potential, 

V  ( r  )  -   qe   /  ™  r 

is  convenient  to  consider  since  Calogero  treats  the  case  with  strength 
parameter,  g  =  -10,  m  =  1.  The  output  produced  by  the  second  program 
is  graphed  in  Figs.  2  and  3. 

The  most  obvious  feature  of  the  graphs  is  the  presence  of  "levels" 
which  arise  from  the  presence  of  possible  bound  states  in  the  potential. 
Calogero  discusses  the  connection  of  this  with  Levinson's  theorem  and 
shows  that  this  structure  in  the  phase  function  can  be  used  to  locate 
bound  states.  This  very  feature  though  poses  a  distinct  problem  to 
the  evaluation  of  the  exact  phase  shift.  Owing  to  the  presence  of  this 
structure,  a  simple  test  for  the  end  of  the  integration  that  depends  on 
the  phase  function  tending  to  a  constant  value  at  infinity  is  useless. 
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Figure  2 
Plot  of  Phase  Function  for  K  =  0.1,  Yukawa  Potential 
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Figure  3 


Plot  of  Phase  Function  for  K  =  1.0,  Yukawa  Potential 
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For  this  reason,  the  integration  must  be  carried  out  past  any  possible 
structure  before  a  cutoff  of  the  above  type  is  allowed.  In  order  that 
the  program  will  do  just  that,  an  integer  counter  is  established  and 
required  to  be  larger  than  a  certain  lower  limit  which  is  designed 
to  carry  the  integration  over  the  structure.  The  procedure  for  doing 
this  is  as  follows;  the  second  form  of  the  program  is  first  run  with 
the  intention  of  finding  the  structure  and  determining  the  value  of  DR 
which  gives  the  best  results.  Having  found  the  structure/  the  counter 
may  be  set  appropriately.  By  use  of  this  artifice,  the  program  is 
enabled  to  reproduce  the  results  given  in  Calcgero's  book  to  within 
a  few  percent.  It  should  be  noted  that  another  way  of  treating  this 
problem  is  to  establish  a  lower  limit  on  the  basis  of  the  value  of  r 
reached.  This  approach  would  be  superior  if  there  were  a  method  for 
predicting  the  position  of  the  last  flat  region  in  n  •  The  counter 
method  was  chosen  only  on  the  basis  of  a  desire  to  know  the  number 
of  iterations  used,  for  computer  timing  purposes. 

A  further  pecularity  with  the  direct  integration  routine  is  the 
sensitivity  to  the  x  increment,  DR.  In  investigations  with  different 

values  of  DR,  it  was  found  that  when 
DR  is  chosen  too  large  the  structure 
at  small  r  tends  to  destroy  the  ac- 
curacy of  the  routine  by  adding 
large  terms  over  some  of  the  flat 
regions  in  the  phase  function.  An 
example  is  for  K  =  0.1  and  dr  =  0.1. 

Fiqure  4  The  result  ^s  *l  =  44.1  in  110  steps 

as  shown  in  the  diagram  (Fig.  4) . 


40     - 
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The  true  result  (in  agreement  with  Calogero's  book,  p.  22)  is 
n  =  6.24  in  38  steps  found  by  taking  dr  =  0.01  or  less.  Now  the 
utility  of  the  second  program  for  preliminary  tests  become  apparent 
since  it  allcws  us  to  establish  the  number  of  bound  states  for  the 
potential  and  the  proper  increment  and  the  lower  counter  limit  in 
just  a  few  test  runs  which  may  be  quickly  run  on  the  computer  time 
sharing  system. 

The  program  also  has  a  very  important  limitation  in  that  the 
accuracy  decreases  sharply  for  a  small  k  due  to  numerical  round  off 
errors  arising  frcm  the  necessity  of  dividing  by  a  small  number.  The 
situation  may  be  improved  by  choosing  a  smaller  increment,  DR,  but 
after  a  certain  point  this  procedure  tends  to  become  self  defeating 
because  decreasing  DR  increases  the  program  running  time.  A  test  run 
on  a  Yukawa  potential  with  g  =  -2,  K  =  0.1,  and  dr  =  0.1  gives  a  close 

answer  to  the  zero  energy  phase  shift,  but  if  the  value  of  k  is  reduced 

-3 
to  K  =  10  ,  the  accuracy  in  the  calculated  value  of  the  phase  shift 

actually  decreases. 
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IV.   RESULTS 

The  final  problem  is  the  application  of  the  programs  to  the 
polarization  potential  discussed  earlier,  eqn.  (12) .  Figures  5 
and  6  detail  the  results  of  the  explicit  calculation  with  the 
second  program.  They  indicate  that  a  value  of  DR  of  0.1  will  carry 
the  program  nicely  through  to  the  asymptotic  region  if  the  lower 
counter  limit  is  set  at  200.  The  lack  of  structure  in  the  graphs 
tends  to  make  this  requirement  useless  though  and  it  becomes  pos- 
sible to  set  the  limit  at  10  iterations  just  to  let  the  program  get 
started  before  the  asymptotic  cutoff  procedure  is  allowed  to  function. 
The  lack  of  structure  just  commented  on  tends  to  indicate  a  lack  of 
bound  states  in  the  potential  as  might  be  suspected  from  the  very 
weak  nature  of  the  function.  Tests  could  be  done  by  imposing  a 

strength  parameter  on  the  potential ,  i.e., 

i 


J(2 


'6 

which  could  then  be  adjusted  as  suited  the  investigator.  By  means 

of  such  an  artifice,  bound  states  might  begin  to  appear  when  the  peak 
of  the  potential  is  high  enough.  In  addition,  increasing  the  strength 
of  the  potential  would  increase  the  magnitude  of  the  phase  shift.  One 
might  expect  that  such  structure  would  appear  for  values  of  r  less  than 
1.0  and  that  the  way  to  locate  such  structure  would  be  to  decrease  the 
value  of  the  increment  DR.  Though  interesting,  such  an  investigation 
does  not  have  any  apparent  physical  use  since  the  only  physical  feature 
of  the  potential  is  the  asymptotic  dependence  at  large  r.  Thus  further 
discussion  of  this  point  has  been  omitted  here. 
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Figure  5 


Plot  of  Phase  Function  for  K  =  0.1,  Polarization  Potential 
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Figure  6 
Plot  of  Phase  Function  for  K  =  1.0,  Polarization  Potential 
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The  most  interesting  feature  of  the  calculation  from  the  stand- 
point of  possible  application  of  the  program  to  real,  physical 
potentials  is  the  necessity  of  carrying  out  the  calculation  to  larger 
values  of  the  radial  coordinate  than  was  the  case  for  the  Yukawa 

potential.  The  reason  for  this  behavior  is  obvious,  the  polarization 

4 
potential  only  drops  off  as  fast  as  1/r  at  infinity.  This  leads  to 

the  observed  difference  in  the  convergence  of  the  scheme  and  probably 
to  the  reason  why  the  only  places  that  the  Phase  method  is  applied  in 
the  literature  are  in  papers  dealing  with  nuclear  type  potentials  which 
may  be  assumed  to  fall  off  very  rapidly.  To  get  some  idea  of  the  magni- 
tude of  this  effect,  the  graph  in  Fig.  2  may  be  examined  in  comparison 
to  the  graph  in  Fig.  5.  The  plot  for  the  Yukawan  potential  has  sta- 
bilized by  r  =  1  while  with  the  polarization  potential,  the  plot  is 
just  stabilizing  at  r  =  25.  These  observations  hold  great  import  for 
any  application  of  the  method  to  atomic  like  potentials  since  the  matter 
of  time  available  for  computation  is  often  very  critical.  Even  on  the 
IBM-360  time  sharing  system,  timing  for  the  program  had  become  critical 
when  increments  less  than  0.1  were  chosen.  In  addition,  output  disk 
space  became  a  problem.  The  problem  of  timing  could  be  expected  to  be 
considerable  worse  on  smaller  computing  systems. 

Having  discussed  the  phase  function  and  its  properties,  use  may 
now  be  made  of  program  1  to  obtain  a  plot  of  n  vs.  energy  as  in  Fig.  7. 
The  most  apparent  feature  of  Fig.  7  is  the  0.78  radian  peak  at  k  =  0.36. 
It  is  a  real  peak  as  opposed  to  an  apparent  peak  that  could  be  produced 
by  a  poor  choice  of  the  increment  at  such  a  low  value  of  k.  An  inter- 
pretation of  this  peak  will  not  be  attempted  at  this  point  but  a  com- 
ment is  appropriate  on  the  behavior  of  n  (K)  on  the  low  energy  portion 
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of  the  peak.  It  seems  to  drop  off  drastically,  possibly  to  zero. 
Below  k  =  0.1,  no  data  taking  was  attempted  due  the  necessity  of 
integrating  too  far  with  too  small  an  increment. 
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Figure  7 


Plot  of  Phase  Shift  vs.  K  for  Polarization  Potential 
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V.   CONCLUSIONS 

The  program  developed  for  the  phase  equation  integration  has  been 
tried  on  a  test  case  and  applied  to  a  second,  semi-physical  case  with 
somewhat  the  expected  results.  This  program  has  in  these  two  cases 
performed  well  as  long  as  the  precautions  mentioned  were  taken.  It 
should  be  noted  that  no  error  analysis  of  the  routine  was  attempted 
due  to  lack  of  time,  however,  the  comparison  of  the  results  on  the 
Yukawa  potential  for  g  =  -10  with  the  results  published  by  Calogero 
in  his  book  [Ref .  7] ,  tends  to  point  to  an  accuracy  that  is  within 
a  few  percent.  Whether  this  result  is  valid  for  any  potential  other 
than  the  Yukawa  potential  and  in  particular  whether  it  is  equally 
true  for  atomic  level  potentials  is  a  question  that  remains  unan- 
swered at  this  point.  It  is  expected  that  simple  numerical  errors 
arising  from  the  Runge-Kutta  integration  scheme  and  the  ordinary 
computer  errors  of  round  off,  etc.  will  decrease  the  accuracy  of  the 
calculation . 
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APPENDIX 

(A.l)  The  Kohn-Hulthen  Method 

The  most  frequently  used  variational  procedure  is  the  Kohn- 
Hulthen  method  [Ref .  15,  16]  which  is  based  on  the  stationary  property 
of  the  functional  A^  *  where 

(Al)  J   =   2m  [  <\>  (  £  -  H  )  +  dr 

tf    J 

and  A  is  some  function  of  the  phase  shift  determined  from  the  normal- 
ization of  V  as  r  —?  &&    .  This  paper  will  employ  tangent  normalization 
throughout ,  i.e., 

(A2)  <fj   >  (l/h)'z       j£<«l0   +   *<t*ify  ntC^r)J 

r—  >^o 

where  fy  is  the  phase  shift  for  the  J?        partial  wave  and   u  and     fy 
are  the  spherical  Bessel  and  Neuman  functions  respectively.     Consider- 
ing the  special  case  of  s-wave  scattering,  Q  =  0,  then   (A2)  becomes 

ip    >       (<■//?)    2  [sm  l<r/«r    +    Von  «<,  coilcr/irrj 

C///)       ccs^  Kr 

With  this  form  for  r  ,  one  now  considers  the  explicit  variation  of  (Al) 

(A3)    51  --  U*kV2  t  k2-  u)4m?-  t-  U(v^<l-0)S'i'<ir 

Since  for  s-wave,  *r  has  no  angular  dependence, 
so  that  the  second  integral  in  (A3)  becomes 
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Integrating  by  parts,  this  gives 


iJr, 


S^i.^-o)\'& 


-f  U^(<!-  uj^d 


oo 


34*  rldj^     i  2^  (^C£"H)^c)i 


Thus  (A3)  becomes 


t    +  l/il 


ao 


0 


The  last  term,  upon  application  of  the  boundary  conditions,  reduces 
to  -  6  0*  A  ^o  /<  )   so  that 

cm,      £tt  +  +«^A)    --   2^v^z-u)^r 

But  the  Schrodinger  equation  is  just 
so  that 

The  stationary  functional,  [A]  =  'Ov^^/jC  ,  is  given  by 

<«>  [M  -"     X    +  I    ■     X  <  l-i  U(£-  HJ^dr 
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The  trial  function,  t    ,   is  chosen  to  have  the  form 

m  I 

where  i  "  vHi/y  ^Jo  ~  ^  K  Ji  ]   and  the  X»   are  some  convenient  basis 

functions,  that  rapidly  approach  zero  as  r-?00  .  The  asymptotic  part 

of  the  wave  function  t  /  is  chosen  so  as  to  shield  the  —  singularity 

in  cos  kr/r  as  r-;  0.  It  is  evident  that  r  has  the  correct  asymptotic 

form  as  -  -  r  /  \  1 

(0      =       (lIu)   'Js'iKr/Kr  -   \k  (itHkr/ictr1-  co^kr  /«r)j 

The  Kohn-Hulthen  principle  then  states  that  the  variation  of   $ ]  with 
respect  to  all  arbitrary  constants ,  i.e., 

(as)  <5[A]  .-   &(\  + 1)  ■•  O 

produces  a  value  of  [A ]  that  is  second  order  accurate  indA  •  The 
explicit  variation  is  done  by  writing  the  functional  in  the  form; 


(A9) 


where 


[A]  ■  I  a  cy  nts   +  a  I  Cj-2|  +  3  +  \ 


h-  = 


2-%    [^CB-HHi     -    fti{ 


B  --     2a  \iCce-M)<J 


<? 
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Variation  of  the  n  constants,  C.  then  gives  the  system  of  linear  equations 
(A10)        Z   M2'j  Cj  *       -   ^i 

Since  the  trial  function  is  linear  in  A  / 
and  since  the  M^i  are  independent  of  A 
This  leaves  the  equations: 

Z  ^  <?/     =    -  Si 

(All)  J 

Solving  for  the  constants/ 
where 

W,-  Zc/'Vfyc/'1  +  2  Jc/°*/" 

Then  a  final  variation  with  respect  to  A  gives,  to  first  order 

\   -   -  (■    W'  ±   B.  ±  /  ) 
and  the  second  order  accurate,  stationary  functional  is 

(A12,      [A]  ■   Vo  .  B0  -    (*/.  +  b,*q* 
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(A.  2)  The  Schwinger  Method 

A  less  frequently  used  but  more  elegant  variational  method  is 
that  due  to  Schwinger  [Ref.  17,  18,  19,  20] .  It  gives  the  stationary 
value  of  the  functional  [\(^otyj  \    directlv.  The  problem  is  to  solve 
the  radial  Schrodinger  Equation. 


(A13)       (f  u ,       \  \  K2  -    2LR±})\  UcCr) 
6\l  [  r2     1 


U(r)uD(r) 


To  do  so,  one  solves  for  the  Green's  Function  Gjj?  (  r  }  f0  \     given  by, 


(A14) 


xT   CiAr  lO    -+    \  Kl-    SLULtH 


dr; 


CiJr  lrc)-   -  <5(r-  r„) 


which  gives 

Cis(rlr0)=         -Krrti    i£(icr  )  n^Cicr*)        r  <  rc 

(A15)  * 

The  integral  equation  for      Uf(r)     is  then 


(A16) 


Ujptr)    -    njj(n)   -     \circ,ClQ(rlrc)UCr0)u£(rc) 


For  large  r,  the  asymptotic  forms  are 

l< r  \8 C \( r )  y    cuz  (tr  -   {(4+  i)n  ) 

Kr   Yij,  (/<rr)  ?     S»m  (wr    -  £  (4i  OuJ 

therefore,  the  asymptotic  form  of  Us(r)  may  be  written  as 


Uo(f) 


?    cov'[k>  -    J  <**  0«J 


(A17) 


.Oo 


+   seller-  i(J?+')il]     <*cv  r0   i,  (*r0)(J(r0}  f<0 

Jo        ^ 
But,  this  is  to  match  to  the  asymptotic  form 

Cr}    —  y    Co*  h>'  -   \  U+l)»  j  -4    z,«  [k  r-  +  Cf*  i)~]  ian  ^ 


Ut 
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and  so  one  sets 


Inserting  the  integral  form  for   Itf(K')  in  terms  of  ^(r)    from  (A16) 


.00 


4a 


d<    <J(^°«tO  \oetr)  t   Grrlr„)  Ufr^W^r.jdfc 


J 


therefore , 


..«5  f       f 


dividing  now  by   "•  a  n  r]ii 


(A18) 


Kco^i,  ^ 


\c\rOCr)ue2Cr)    t     (cir    \  cJf*  U(r  )utic  )  GutHrJOWwt,) 

Jo  Jo  iu 


iU IV 


•&0 


cJr   r    U(r  )  us(r)u({r) 


For  s-waves  this  becomes 


(A19)    Kcdi 


lo 


\c\t  L(r)v*<r>   -    \c\r   \c}r0U(r)»cr)£i(r\r0)U(r0)o(re,') 

)i~-  _  .      ,^  U         Jo  ._ 


-    \dr    UCr)  uLr  )  <,<m  if  r 


where    G?  ^  I  rc )    is  now : 


G^Cr  I  f^     -        ~    Sink  r  cos  k  r0 


r<   r. 


K 
In  order  to  show  that  (A19)  is  in  fact  stationary  to  second  order  for 

the  exact  wave  function,  consider 


(A20) 


cSUco-f  >(]   -[Uchj(,tdu)j  -  [lUoi^U.)] 
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Expanding  and  collecting  terms,  we  have 


Q<  0(f) vlf)  \»(ti  -  \dt0Hnu)  U(f0)u(O\  . 

jo l if +   oCo  ) 

combining  the  terms  the  numerator  for  the  expression  will  be 


-|  2  Wr  t-Kr^  2oCr) 


r 

u(r)  - 

Jo 


dr0  Ci(r\r0)0(To}v<tc) 

4  2[jUu(r)^Ji^»  ^r  j"^|cirOCrJuCr)5^^Jj-/^cyy 


However 


(A21) 


Oo(r0)h 


rW 


kco+r)^  |r  drCXf^^)s^^ri 


'A  U" 


but  frcm  (A17)  this  is 


.Qo 


-o  _ 

After  factoring  out    -  fdr  U(r)u(rM»nkrl  ,   the  numerator  becomes 
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r<X> 


However,  combining  terms,  this  is, 

+  as*) 

But  from  (A16) ,  for  u(r)  the  correct  wave  function, 

ocr  )  -  %m\cr   -  Wfo  0](r\^)uCro)u(ro)   H  O 

Thus,  neglecting  terms  of  order  <5  ,  the  numerator  is  zero  and  hence 

Therefore,  we  see  that  (neglecting  second  order  variations)  the  expres- 
sion is  stationary  when  u  is  the  correct  solution  of  the  Schrcdinger 
Equation. 

Now  going  back  to  equation  (A19)  and  inverting 


(A24) 


l*k«?J  = 


-  le „      „. .4.. 


'oo       roc 


dt  Oir)uhr)  -  \6r  \c\r0{jtf)oir)^CrlY0)LXrc)otr^ 

J  0  Jo         Jo 

Letting  o(r )  ~  c^i  %i       then 

(l^v,)5  _          Zcj  c-,  Mj 


(KB)[\L^n%\:    = 


where 


and 


/-oo 

M"  5     \clr  iKr)  j|^rjs/**r 
Jo 

C7lj  =     U    (dV„  Ut«-)^^r)^ir;r0)0^Cr0)-C7ji 
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The  problem  remaining  is  to  solve  for  Y K "»<* ^ %J  and  the  C  5  .  To 
do  so,  define 

(A26)    Q  =    1  CyCj  ([K+»*fy]Mfj  -  A^  )  E  O 

but  Q  is  obviously  stationary  and  so, 
(A27)       3Q   =   2lci([^a^fl]^'-A{-/)-'0 

In  order  to  have  a  solution  for  the  C^  /  the  determinant  of  the  matrix 
of  coefficients  must  be  zero,  i.e., 


(A28) 


c!e^     ([IC +*ify]rtij-AJij  )  -  C 


However,  it  is  a  simple  matter  to  show  that  the  matrix  N  is  of  rank 
one,  i.e.,     there  are  no  non-zero  sub-determinants  of  the  matrix  N  of 
order  2,  and  further,  it  can  be  shown  that  such  a  system  has  only  one 
nontrivial  solution  for      [_  K  •  #  A  f]0  [  , 


(A29) 


/  Mi.     "ti   >  - 
Mi,    Mzz    •  ' 


[Kh*$  " 


i 

Mm 
/ 

M,2 
A/12 

A 

•     /      f 


Thus  one  finds  the  solution  for  J  K  "f  0  w  J9  /  without  recourse  to  a 
solution  for  the  C3  .  This  result  is  shown  in  Kohn's  Paper  [Ref .  15] . 
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